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Accumulating evidence for nonstandard leptonic decays of D s mesons 
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The measured rate for ^£ + v decays, where £ is a muon or tau, is larger than the standard 
model prediction, which relies on lattice QCD, at the 3.8a level. We discuss how robust the the- 
oretical prediction is, and we show that the discrepancy with experiment may be explained by a 
charged Higgs boson or a leptoquark. 

PACS numbers: 13.20.Fc,12.60-i,14.80.-j 



Introduction. — The pattern of flavor and CP violation 
of the standard model has been established by a wide 
range of experiments. This agreement, however, leaves 
room for new flavor effects to show up as calculations 
and measurements improve. Intriguingly, decays of the 
D s meson (the lightest cs state) could be more sensitive 
to new physics than any other process explored so far. It 
suffices that a new particle couples predominantly to lep- 
tons and up-type quarks, but not to the first generation. 

In this Letter we examine the leptonic decays of the D s . 
Recently, the calculation of the relevant QCD matrix 
element has improved significantly, and more accurate 
measurements of the rate have been made. The average 
of the experimental results disagrees with the standard 
model by almost four standard deviations. We discuss 
the evidence, and propose that a nonstandard amplitude 
interferes with the standard IL-mcdiatcd amplitude. We 
show that the tree-level exchange of a spin-0 particle with 
mass of order 1 TeV may account for the discrepancy. 

Leptonic D s decays. — The D s —>tv branching fraction, 
where £ is a charged lepton of mass m#, is given in the 
standard model by 



B(D s ^tv) 



m Ds 
8tt 



TD.fb. \G F V: s m e \ z 1 



(1) 



Here mo, and td 3 are the mass and lifetime of the D s , 
Gf is the Fermi constant, and V cs is a Cabibbo- 
Kobayashi-Maskawa (CKM) clement. The decay con- 
stant fo g is defined by 



(0| s-f^c \D s (p)) = if D „Pn, 



(2) 



where p^ is the 4-momentum of the D s meson. Al- 
though the electroweak transition proceeds at the tree 
level, Df — > W + -^t + v e , its rate is suppressed. The he- 
licity of the lepton must flip, leading to the factor mi in 
the amplitude. For the muon, this helicity suppression 
(me/mD e ) 2 is 2.8 x 10~ 3 . The r mass is only 10% smaller 
than the D s mass (1.969 GeV), so there is no significant 
helicity suppression, but the phase space suppression [the 
last factor in Eq. H}] is 3 - 4 x 10 ~ 2 - 



TABLE I: Experimental values of fr> a ■ Our averages treat 
systematic uncertainties as uncorrelated and omit the PDG 
entry [l[, which is an average of earlier experiments. 



final state 


reference 


So. (MeV) 


£v 


PDG [1] 


294 ± 27 


1"' 


BaBar [3] 


283 ± 17 ± 16 


/"' 


CLEO [4] 


264 ±15± 7 


flV 


Belle [5] 


275 ± 16 ± 12 


TV (T^ITV) 


CLEO [4] 


310±25± 8 


tv (r—^evv) 


CLEO [6] 


273 ±16± 8 


flV 


our average 


273 ± 11 


TV 


our average 


285 ± 15 



We have collected in Table Q] all precise experimental 
measurements of B(D* — > tv), which are usually quoted 
in terms of fr> s [H, 0|. Combining the error bars in 
quadrature, our average of tv and [iv final states is 

(/oJex P t= 2 77±9MeV. (3) 

The most accurate calculation from lattice QCD is 0] 

lfrJ Q CD = 241±3MeV, (4) 

where statistical and systematic uncertainties arc com- 
bined in the fitting methods. The only other modern 
lattice-QCD calculation agrees, 249 ± 3 ± 16 MeV 0], 
but its quoted error is five times larger and would not 
influence a weighted average with Eq. (|4]). The discrep- 
ancy between Eqs. © and g|) is 15% and 3.8cr. Tabic U 
also shows averages for each mode separately: for tv (fiv) 
alone, the discrepancy is 18% and 2.9a (13% and 2.7a). 

If the BaBar result is omitted from the average, as in 
Ref. 0, then the discrepancy is 3.4er. On the other hand, 
if the earlier measurements [1( as well as the BaBar result 
are included, we find a 4.1er discrepancy. 

Experiments. — CLEO 0, Q produces D s pairs near 
threshold, where the multiplicity is low. Their method 
reconstructs one D y s 1 and then counts how often the 
opposite-side D s decays lcptonically. When the charged 
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lepton is a muon, the neutrino is "detected" by requir- 
ing the missing mass-squared to peak at zero. When the 
charged lepton is a r, the identification is made through 
the subsequent decays r — ► evv and t — ► irv. BaBar [3[ 
observes D s coming from the decay D* ^D s 7, produced 
well above threshold. They compare the relative num- 
ber of subsequent D s — ► ji + v and D s — * cjm, and then 
use their own measurement of B(D s —Kf>Tr) to determine 
B{D s ^iv). Belle 0] also observes D s via D* ^ D s j, 
but the whole event is reconstructed, using a Monte Carlo 
technique. In summary, all these measurements have cen- 
tral values and error bars that are straightforward to in- 
terpret, and to combine to obtain Eq. (|3j). 

The measured branching fraction and Eq. ([1]) yield 
\Vcs\fD s - Three-generation CKM unitarity is assumed, 
either taking \V CS \ from a global fit to flavor physics 
or setting \V CS \ = \V U d\- The difference is numerically ir- 
relevant. Relaxing the assumption cannot lead to agree- 
ment between theory and experiment because unitarity, 
even for more than three generations, requires |y cs | < 1, 
whereas the discrepancy would require |V CS | » 1.1. 

Radiative corrections. — The measurements are not, 
strictly speaking, for D s — > tv alone, because some pho- 
tons are always radiated. The radiative corrections have 
been studied, focusing on effects that could overcome the 
helicity suppression jg. fl~oj] - 

For D s ~^t + v there is no sizable helicity suppression. 
In the rest frame of the D s , the r acquires only 9.3 MeV 
of kinetic energy, so it cannot radiate much. Explicit 
calculation [§] shows that the radiative corrections are 
too small to account for the discrepancy fill ]. 

For D s — > fi + i> radiative corrections could play a role 
due to processes of the form D s — > 7D* where D* 

is a (virtual) vector or axial- vector meson. The transition 
D* — > n + v is not helicity-suppressed, so the factor a for 
radiation is compensated by a relative factor rn 2 Ds /m 2 ^ for 
omitting helicity suppression. Using Eq. (12) of Ref. § 
and imposing the CLEO Q cut E 7 > 300 MeV, wc find 
that the radiative rate is around 1% and, hence, insuffi- 
cient to explain the discrepancy. 

Lattice QCD — There are many lattice-QCD calcula- 
tions for fr> s in the literature, but only Refs. 0, 0] in- 
clude 2+1 flavors of sea quarks, which is necessary to find 
agreement for many "gold-plated" quantities, namely 
those for which errors are easiest to control 12j. Both 



calculations start with lattice gauge fields generated by 
the MILC Collaboration [3] , which employ "rooted stag- 
gered fermions" for the sea quarks. At finite lattice spac- 
ing this approach has small violations of unitarity and lo- 
cality. Theoretical and numerical evidence suggests that 
these vanish in the continuum limit, such that QCD is ob- 
tained, with the undesirable features controlled with chi- 
ral perturbation theory. The strengths and weaknesses 
of this approach have been reviewed in detail (l3 |. 

Reference [7j reports an error five times smaller than 



that of Ref. [8| for several reasons. The largest uncer- 
tainties in Ref. Q come from a power-counting estimate 
of the discretization error for the charm quark, and from 
uncertainties in the chiral extrapolation. Reference 
employs a different discretization for the charm quark, 
which allows a controlled extrapolation to the continuum 
limit. Thus, the discretization error here is driven by the 
underlying numerical data. 

The action for the charm quark in Ref. [7|, called 
HISQ [15j, is the same as that used for the light valence 
quarks. As a result the statistical errors are smaller than 
those of the heavy-quark method used in Ref. 0], and 
the axial current automatically has the physical normal- 
ization. The suitability of HISQ for charm is one of its 
design features, it has been tested via the charmonium 
spectrum [l5| , and the computed D and D s masses agree 
with experiment. The D + decay constant f D + also agrees 
with experiment, at lcr. 

Another feature of Ref. [7| is the way the lattice- 
spacing and sea-quark mass dependence is fitted. Full 
details are not yet published, but it is noteworthy that 
the same analysis yields f n and Jk in agreement with 
experiment [l| and earlier, equally precise, lattice-QCD 
calculations [16[ . The D s meson is simpler than the pion 
or kaon for lattice QCD, because none of the valence 
quarks is light, so fn s is easier to determine than / w . We 
find that simple extrapolations lead to the same central 
values for both to_d 3 and fjj a . 

The error bar in Eq. (j4]) is smaller than that in Eq. |3]). 
Therefore, it is the combined experimental error that pro- 
vides the yardstick for the deviation. To illustrate, if 
the lattice-QCD error bar were doubled, the discrepancy 
becomes 2.7a, 2.5a, and 3.3er for t, /i, and combined. 
Hence, even if additional sources of uncertainty are un- 
covered, evidence for a deviation may well remain. 

Nonstandard effective interactions. — Although the ex- 
periments quote the final states as /x+z/^ and t + v t (and 
their charge conjugates), the flavor of the neutrino is not 
detected. Nonstandard physics could lead to any neu- 
trino flavor, even a sterile neutrino. However, given the 
large effect that needs to be explained, we shall restrict 
our attention to amplitudes that could interfere with the 
standard model, which fixes the neutrino flavor. Lorentz- 
invariant new physics may contribute to D s — * ivi only 
through the following effective Lagrangian: 



C 



^(s l5 c)(D L e R )+U.c, (5) 



where C\ and Cp are complex dimensionless parame- 
ters, M is the mass of some particle whose exchange in- 
duces the 4-fermion operators ([5]), and the c,s,£ fields are 
taken in the mass-eigenstate basis. The hadronic matrix 
element required for the decay induced by (s~J5c)(vlIr) 
is related to the one of Eq. ([5]) by partial conservation 
of the axial current: (ra c + m s )(0\ sij§c \D S ) = fDs'rnjj ■ 
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The branching fraction in the presence of the operators 
([5]) is given by Eq. ([!]) with GpV* s me replaced by 



G F V* s me 



1 



V2M 2 



C A me, 



Cp m 



(6) 



with no helicity suppression in the last term. 

The imaginary part of V cs is negligible (in the standard 
CKM parametrization so constructive interference, 
which would increase B(D s —>£v), requires the real part 
of C A or Cp to be nonzero and positive. Assuming only 
one nonzero coefficient, the amplitude for t + v t (/i + v^) 
could be increased by 12% (8.4%) only if 



M 



(ReC A y/ 2 
M 

(ReCp) 1 / 2 



< 



< 



710 GeV for £ = r 
850 GeV for £ = n 

920 GeV for £ = r 
4500 GeV for £ = fj, 



(7) 



(8) 



thereby reducing the discrepancy to la in each case. 
These bounds are a key new result of this Letter, be- 
cause they constrain any model of new physics. 

The effective interaction also contributes to the 
semileptonic decays D — > K[i + v. This proceeds through 
two amplitudes, corresponding to angular momentum 
J = 1 or for the lepton pair. For J = 1, the standard- 
model amplitude and that from C A are not helicity sup- 
pressed, while that from Cp is. For J — 0, the pattern of 
helicity suppression is as for the leptonic decay. Hence, 
only the J = 1 part of the rate will be visible, and as 
the accuracy of the lattice-QCD calculations improves, 
the comparison with experiment will help decide which 
interactions are responsible for the effect in D s — » £v. 
The current status favors Cp ^ rather than C A ^ 0, 



because the lattice-QCD 
agrees with experiment [If 



irediction for D 
L albeit at the ~ 



— > K/jv 
7% level. 



17] 



New particles. — There are three choices for the electric 
charge of a boson that can mediate the four-fermion op- 
erators ([5]): +1, +2/3, —1/3, corresponding to the three 
diagrams shown in Fig. [TJ The exchanged boson (taken 
to be emitted from the vertex where c is absorbed) is 
a color singlet if the electric charge is +1, and a color 
triplet if the electric charge is +2/3 or —1/3. We shall 
consider only the cases where the new boson has spin 
or 1, and its interactions are renormalizable. 

A new vector boson, W' , of electric charge +1 would 
contribute only to C A . Such a boson must be associated 



(+1) 




£+ 



t (+2/3) t (-1/3) 



«. £+ '' 
FIG. 1: Four-fermion operators induced by boson exchange. 



with a new gauge symmetry, which makes it difficult to 
allow large couplings to left-handed leptons. One pos- 
sibility is that W and W mix, but the constraint from 
electroweak data on mixing (< 10 -2 ) is too strong to 
allow noticeable deviations in D s decays. Another pos- 
sibility is that some new vector-like fermions transform 
under the new gauge symmetry and mix with the left- 
handed leptons. Such mixing is also tightly constrained, 
especially by the nonobservation of vector-like fermions 
at LEP and the Tevatron. Overall, a W is inconsistent 
with Eq. ([7]) , barring perhaps some finely-tuned elaborate 
model (e.g., with large W-W mixing whose electroweak 
effects are cancelled by other particles). 

A spin-0 particle of charge +1, H + , appears in models 
with two or more Higgs doublets. Its interactions, in the 
mass eigenstate basis for charged fermions, include 



H + {vcCrs l + y s c L s R + yivit) + H.c, 



(9) 



where y c ,y 8 ,yi are complex Yukawa couplings. The ex- 
change of H + induces C A = and 



(10) 



taking M equal to the H + mass. If H + is the charged 
Higgs boson present in the Type-II two-Higgs-doublet 
model, then y c Jy s = m c / (m s tan 2 /?) so that Cp can have 
either sign [l9(, but the Yukawa couplings are too small 
to be compatible with Eq. ijHJ). Other models lead to large 
constructive interference. For example, a two-Higgs- 
doublet model where one doublet gives the c, u (but not 
d, s, b, or i) and lepton masses, and has a vacuum ex- 
pectation value of ~ 2 GeV, yields \y s \ <§; y T , y* ~ 0(1). 
Thus, Cp > and the limits © are satisfied for M < 500 
GeV. Furthermore, such a model explains why the devia- 
tions in tv and /iv are comparable. It is encouraging that 
this two-Higgs-doublet model does not induce tree-level 
flavor-changing neutral currents, and the off-diagonal 
couplings of H + are CKM suppressed. Given that this 
model has not been previously studied, its 1-loop con- 
tributions to flavor-changing processes (such as b — > sj) 
need to be computed before deciding whether some fine 
tuning is required to evade experimental bounds. 

The charge —1/3 and +2/3 exchanges correspond to 
leptoquarks. A scalar charge +2/3 exchange arises for 
the (3, 2, +7/6) set of SU{3) c xSU{2) w xU(l) Y charges. 
This leptoquark appears, for example, in a new theory 
of quark and lepton masses [icj. Let r — (r u , r^) be the 
doublet leptoquark, where ra is its charge +2/3 compo- 
nent. The interaction terms relevant here, written in the 
same basis as ((SJ, are \ c iTdCRV l L + X' se rdSL£R- The 
exchange gives C A = and C P = — A* £ A^/4. Since the 
leptoquark couplings can have any phase, the new ampli- 
tude can interfere constructively. Still, various flavor pro- 
cesses constrain the couplings of r. Even if its couplings 
to first-generation fermions were negligible, the lepton- 
flavor violating decays r — > /xss, where ss hadronizes to 
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77, n' , </> or KK, set a lower limit on M 2 /|Ag T A' |, which 
is hard to reconcile with Eq. ([8]). One way out would be 
a model with two r leptoquarks, with one coupling to r 
and the other one to ji. The constraint from t — * /j,ss 
similarly disfavors spin-1 leptoquarks of charge +2/3. 

A scalar leptoquark of charge —1/3 (also discussed in 
[20| ) arises in the case of two sets of SU(3) C x SU (2)w X 
U{1)y charges: (3, 1, —1/3) or (3, 3, —1/3). Let us denote 
the former by d. Its Yukawa couplings are given by 



Fermilab is operated by Fermi Research Alliance, LLC, 
under US DoE Contract DE-AC02-07CH11359. 



d [m (c L £ c L - s L v e L c ) + 4 c R £ c R ] + H.c. 



(11) 



where Kg and n' e are complex parameters. These inter- 
actions are present, for example, in i?-parity violating 
supersymmetric models (their effect on D s —>■ e + v has 
been analyzed in Ref. [2lJ). The d exchange, as in the 
last diagram of Fig. [IJ gives (for M equal to the d mass) 



C\ = 



(12) 



For \K^/nt\ <C mim c /m 2 D , the interference is automati- 
cally constructive [see Eq. (jBJ)], and the resulting devia- 
tions in tv and \xv are approximately equal if ss |/c t |. 
Moreover, there are no severe constraints from other pro- 
cesses on the couplings ki and n' e with £ = t or /j. The d 
couplings to the electron can be forbidden by a symmetry, 
and the ones to first-generation quarks could be small. 

The (3, 3, —1/3) scalar leptoquark includes an SU(2)w 
component of charge —4/3 which mediates r— >/jss. The 
vector leptoquark of charge —1/3 has the same problem. 

Conclusions. — We have argued that the 3.8<r discrep- 
ancy between the standard model and the combined ex- 
perimental measurements of D s — > Iv appears so far to 
be robust, and thus it is worth interpreting it in terms 
of new physics. The upper bounds ([7]) and ([8]) on the 
scale of four-fermion operators are low enough to allow 
exploration of the underlying physics at the LHC. 

A d scalar leptoquark of charge —1/3 may solve the D s 
puzzle without running into conflict with any other mea- 
surements. At the LHC, the d can be strongly produced 
in pairs, and the final states would be £ + £~ jj, where £ is 
a t or a ji, and j is a c-jet. Given that there are two £j 
pairs, each of them forming a resonance at the d mass, 
the backgrounds can be kept under control. The current 
limits on the d mass from similar searches at the Tevatron 
are around 200 GeV Q. 

An alternative explanation is provided by an H + ex- 
change in a (new) model where a Higgs doublet gives 
masses to the charged leptons and c and u quarks, and 
a second Higgs doublet gives masses to the down-type 
and top quarks. Both the leptoquark and charged Higgs 
solutions lead naturally to comparable increases in the 
branching fractions for D s ^t + v and D s — >fi + v, as sug- 
gested by the data. 
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